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Abstract 

We consider a dielectric medium with an ultraviolet behavior as it 
follows from the Drude model. Compared with dilute models, this has 
the advantage that, for large frequencies, two different media behave the 
same way. As a result one expects the Casimir energy to contain less 
divergencies than for the dilute media approximation. We show that the 
Casimir energy of a spherical dielectric ball contains just one divergent 
term, a volume one, which can be renormalized by introducing a contact 
term analogous to the volume energy counterterm needed in bag models. 
PACS: 12.20.Ds, 03.70.+k, 77.22. Ch 



Due to recent progress in experimental techniques, it is now possible to 
measure Casimir forces between macroscopic testbodies very accurately ||. 
Even quantum mechanical actuation of microelectromechanical systems by the 
Casimir force is under discussion ||]. For these reasons a thorough theoretical 
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understanding of the Casimir effect witli realistic media is desirable. There 
are essentially two ways to calculate the Casimir energy. One possibility is to 
sum up retarded van der Waals forces between individual molecules |^ . The 
second way makes use of quantum field theory in the background of a dielectric 
medium (see e.g. [||, ^, |9[ |l^, |l^, |l^). In this second approach one 
tries to recover retarded van der Waals forces by calculating the vacuum energy 
for the electro-magnetic field in a dielectric background. The relation between 
these two approaches is not well established. Only for a dilute ball, up to the 
second order of a perturbative expansion, the methods have been shown to yield 
the same answers In all other cases of dielectric media the presence of 

divergencies forbids to give a physical interpretation to the finite parts obtained. 
Even for a dilute medium, once the full dependence of the pole on the dielectric 
constants is considered, the result is not even expressible in terms of known 
special functions |16|. So the counterterms needed to renormalize the divergent 
contributions is extremely complicated and the interpretation of the classical 
model associated with the couterterms is completely unclear. 

One possible reason for this briefly described difficulty might be the fact 
that the models of dielectric media treated in the quantum field theory frame- 
work mostly do not fulfill a realistic frequency dependent dispersion relation. A 
direct consequence of an unreasonable ultraviolet behavior is, for example, the 
appearance of divergencies which make it impossible to extract, in a physically 
reasonable way, a finite value for the Casimir energy. 

For those reasons, it seems natural to analyze the divergencies of the Casimir 
effect in a model with a dispersion relation that shows realistic features, at least 
in the ultraviolet range. As we will see, once this behavior is assumed, the 
pole structure resulting in the ^-function regularization scheme is in fact very 
simple. Namely just one pole exists and it is a term completely analogous to the 
one renormalized via the bag constant in the bag model. So, introducing only 
one counterterm, namely a volume energy, the Casimir energy can be rendered 
finite. Apparently there is no intrinsic way within the field theory framework, to 
fix the finite part of this counterterm. In principle, as is the case in bag models, 
it has to be fixed by experiments. However, given the simple structure of the 
renormalization needed, we feel this is an important step in order to understand 
the relation between different approaches employed and to extract finite Casimir 
energies in realistic media. 

So, let us start with a description of the simple model we are going to ana- 
lyze. We will consider a nonmagnetic (/i = /io) dielectric ball of radius a, with 
a frequency dependent permittivity given by e(a;) = eo(l — ^), immersed in 
another medium of permittivity eg. This frequency dependent permittivity is 
the so-called plasma model, which follows from the Drude model in the high fre- 
quency approximation. The phenomenological parameter f2 is usually referred 
to as the effective plasma frequency. Given that we are going to concentrate 
only on the pole structure of the Casimir energy for the described setting, this 
high frequency model is sufficient for our purposes. Formally, the Casimir en- 
ergy of this configuration is defined by summing over the vacuum energies of 
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each mode of the electromagnetic field, 



Cas 



(a) =^ 



Tic 
2a 



(1) 



where the dimensionless quantities z„ — aujn/c are the eigenvalues associated 
with a radius one ball. As it stands, the sum in Eq. (|l|) is divergent, and 
a regularization procedure must be adopted. To properly define the Casimir 
energy we will use the zeta function regularization. In this scheme, Ecas{o) is 
defined as the analytic continuation of the function 
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he 

2a 



(2) 



where ({s) is defined through the eigenfrequencies of the electromagnetic field 
by means of the series 



(3) 



which is absolutely and uniformly convergent for 5R(s) sufficiently large. As is 
well known, usually the zeta function has a pole at s — —1, in which case a 
suitable rcnormalization procedure has to be employed. 

In Eq. (^), the eigenfrequencies w„ must be determined by solving the field 
equation 



AE +ne— 
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(and similarly for if), subject to matching conditions appropriate to the inter- 
phase between two dielectric media. 



Ee, 



Ed., 



1 

Ail 



1 



(5) 



One can consider the transversal electric (TE) modes, for which the electric 
field has the form 

'Si,„r=Mr)tn„^{0,<P), (6) 

and separately the transversal magnetic modes (TM), with the magnetic field 
given by 

'^L„,=gi{r)tYi^„,{e,4>). (7) 



Here, as usual. 

The imposition of the conditions in Eq. (|^) leads to 

Mr = a+)=Mr = a-), 



(8) 



(9) 
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for the TE modes, and 

gi{r = a+) = gi{r = a"), 



— dr [rgi{r)] 
eo 



dr[rgi{r)] 



(10) 



for the TM modes. 

In order to have a discrete spectrum, one might encfose the system into a 
large conducting sphere of radius R, sending i? ^ cx3 at a suitable intermediate 
step. This implies also the following boundary condition at r = i? for fi{r) and 
giir)- 

fiir)UR = 0, dr{rgiir))l^j, = 0. (11) 

Alternatively it is possible to use a formulation in terms of the Jost function 
of the corresponding scattering problem, where a subtraction of the Minkowski 
space contribution is performed at the beginning and the infinite volume limit 
taken impHcitly fl^ . 

It is straightforward to see that the cigcnfrcquencies for TE modes are de- 
termined by the zeroes of the function (we put v = I + 1/2) 



where 



<Jl{zi) {yu{za)Mz2) " Mza)y.{z2)} ., 



Ju{w) = W ji{w) = ^1 —Juiw) 



(12) 



yt,{w) = w yi{w) = J —Y^{w) 



(13) 

are the Riccati - Bessel functions, being z — a{uj/c), z\ = zn{z), Z2 = ztiq, 
zq = zRriQ/a, ^ = n(z)/no, no = ^/e^, and n{z) = ^/e{uJ). 

Similarly, for the TM modes the eigenfrequencies are determined by the 
zeroes of 

AP^(z) = Mz,) {yl{z2)X{zo) - X{-Z2)y'Jzo)} - 

1 (14) 

-^JU^i){yAz2)JU^o)-JA^2)yl{zo)}. 

We will suppose that these functions have only real and simple zeroes in the 
open right half plane. 

In the following we will consider explicitly only the TE modes, since the 
treatment of the TM modes is entirely similar. Due to the spherical symmetry 
of the problem, the ^-function has the general appearance 

oo 

i/=3/2 n 
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where v labels the angular momentum and, for a given the index n labels the 
zeroes of A^-^(z) in Eq. (p^). We need to construct the analytic continuation 
of this function to s w — 1. As we will see, C(s) has a simple pole at s = — 1 
and Ete{o-), defined as in Eq. (|^), contains a divergent term which depends on 
the particular dispersion relation adopted for e(a;). The method employed for 
the analytic continuation has been explained in great detail, e.g., in 
|l9| , pO| , and we will simply follow this procedure. However, let us mention that 
the idea of applying this method in the specific context of dielectric media goes 



back at least to 21 



For 3?(s) large enough, we can represent the ^-function as an integral in the 
complex z-plane employing the Cauchy's theorem. For the TE modes we have 



°° 1 f \TE'( \ 

a.):=Ea = ^£^-44f^^^ (16) 

where the curve C encloses counterclockwise all the positive zeros of A^^(z). 
This curve can be deformed into a straight vertical line, crossing the horizontal 
axis at 3?(z) — x, where x is any value satisfying < x < z^^i, Zi,,\ being the 
first positive zero of A^^ (z) . Obviously the integral in eq. ([l6| ) does not depend 
on the particular value of x in this range. 

Indeed, expressing the integrand in terms of the modified Bessel functions 
and taking into account their asymptotic behavior for large arguments, it is 
easily seen that the integral 

converges absolutely and uniformly to an analytic function in the open half-line 
s > 1, which can be meromorphically extended to the whole complex s-plane. 
Expression ( |l7|) can also be written as 

^^^^) = —^^\- - ^l^y-'^ afm^^'J' ^ ^ 

where the prime means derivative with respect to the argument, and we have 
made use of the properties of the Bessel functions of complex argument. 

Changing the integration variable to t = z{y ^ t), with z = > 0, we 
finally get 

1 ^ / + 1) r''\-s A^^ M)) 



C,{s)^~-ny~^e 2^ 'J t-^^ ^—^^dtj. (19) 

Notice that the right hand side of Eq. ( [T9| ) does not depend on z for z(> 0) 
small enough. 
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In order to construct the analytic extension of the expression in Eq. (g) to 
s ~ — 1, we subtract and add the first few terms of the asymptotic expansion of 
the integrand in (^9|) (obtained from the uniform asymptotic Debye expansion 
for the modified Bessel functions appearing in A^-^ [ivt)). In fact, in order to 
isolate the singularities of the Casimir energy it is sufficient to retain in this 
expansion terms up to the order v^'^: 

*^|iM).D-W + <3(.-). (20) 

with 

(21) 

where D^^(i), k — 1, —3, are algebraic functions of t given in the Appendix. 
Notice that we have also discarded terms which are exponentially vanishing for 

i? ^ CXI. 



So, we must consider the series 

TT 




Di''{t) dt ) - 



i/=3/2 



(22) 

The second term in the right hand side of ( P4 ) converges for s > —2 by con- 
struction. Therefore, we can put s = — 1 inside the sum and the integral, and 
numerically evaluate this contribution when necessary. 

In order to investigate the ultraviolet divergencies in Ecas it is sufficient to 
consider the first term in the right hand side of (p3)n. 



^ It should be stressed that the real part in the argument of the series in the first term in 
the right hand side of Eq. (02) is in fact independent of z: Taking into account the analyticity 
of the integrand, one can investigate the 2-dependence by studying the integral 

1 

t-" D'^'^it) dt, (23) 

■tz 

which can be exactly solved in terms of hypergeometric functions. It is straightforward to 
verify that the 2— dependent par t, is imaginary for all s > 1, and therefore is dropped out 
when taking the real part in Eq. (E2). This feature will be useful in what follows. 
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For s > 1 we can study each term in D^^{t) separately, and evaluate the 
expressions 

J2 ^-(-'=-1) / D^^%it) dt , (24) 

1^=3/2 [ ^ -^^^ J 

with the index /c = 1, — 3 corresponding to the order in the Debye expansion. 

Since the real part of the expression between brackets in ( |2^ ) is independent 
of V (it is independent of z = x/v - see footnote ^), the sum over v can be 
performed by means of the Hurwitz zeta function, 

oo 

^ ^-{s-k-D ^ (-^(^^ - /c - 1, 1/2) - 2'~'''\ (25) 

1^=3/2 

Only for k = —3 one gets a singularity at s = — 1 given by 

Ch{s + 2, 1/2)|^^_, = ^ + (7 + ln(4)) + 0{s + 1). (26) 

For the other values of fc a regular extension to s = — 1 is obtained. 

Therefore, it is only for k = —3 that we need to calculate both finite and 
singular parts of the analytic extension of the integral in Eq. ( p^ around s = — 1, 
while for the other values of k only the singular terms are needed. 

All these terms can be worked out exactly to give for the divergent part of 
the contribution of the TE modes to the Casimir energy (in the limit R oo), 

ETEia) = ^ . - . . + Ois + lA . (27) 



a { 24 7rc4 (1 + s) 47rc2(l + s) 

Notice that the first term in the right hand side corresponds to a volume con- 
tribution, while the second one is a curvature contribution. In agreement with 
[ p2| , no surface contribution appears. 
For the TM modes one gets similarly 

ETMia) = ^ , + -^^1^ + Ois + lA . (28) 

a { 24 7rc4(l + s) 47rc2(l + s) ' j 

Finally, adding up the contributions coming from the TE and TM modes, 
we get for the singular piece of the Casimir energy of the dielectric ball just a 
volume contribution, 

since the divergent curvature terms have canceled out between TE and TM 
modes. Replacing the volume of the ball, V — 47ra'^/3, we find 



7 



Thus, we have found that the pole structure of the Casimir energy in the 
framework of the ^-function regularization scheme is very simple, once a realis- 
tic behavior of the dielectric medium at high frequencies is assumed. The pure 
volume divergence might be seen to represent a contribution to the mass density 
of the material and it is to be absorbed into a phenomenological countert- 
erm, much in the way it is done in bag models. Neither a surface tension nor a 
curvature counterterm is needed in this model, since divergencies nicely cancel 
out between TE and TM modes. 

Given this simple pole structure it makes sense to analyze further finite parts 
of the Casimir energy in a quantum field theory context for realistic media. As 
a first step one might assume the Drude model for the medium and analyze how 
the Casimir energy depends on e.g. the plasma and the relaxation frequency. 
More ambitiously, the dielectric constants might be obtained from tabulated 
refractive indices using the Kramer-Kronig relation [23| . This is under consid- 
eration and will allow a detailed analysis of the Casimir energy for a spherically 
symmetric situation with realistic media. 
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A Debye expansion for the TE and TM modes 

The functions Di^{t), k = 1, —3, coefficients of the expansion D'^^{t) (see 
Eqs. ( pO| ) and (|2l|)), which are obtained from the uniform asymptotic Debye 
expansion for the modified Bessel functions appearing in {ii^t), are given 

by 



4")(0^(2^+ ^"°^f ' , (32) 
2a2 (l + nWy 
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a 



4^2 Ano^RH^Z^ no^RH^Z^ 



(34) 

and 

(_3) _ no^ R'^t [64 a'^ -560 a* no'^ +i56 no^ R'^t* -25 no'^ R'^t'^) 

128 I i_|_ "Q"« * \ 

^ ° ^ (35) 

(16 Z^+56no^ Z^+6no*t^ Z^+Sno'**'' Z^)-2no* (t^-SSt* Z^)) 

16*5 (l+no2t2)5 ' 

Similarly, for the TM modes the functions D^l[{t), fc = 1, — 3, appearing 
in the expansion (t) , are given by 

(t) - V ^ (36) 

n(o) (+\ - zl 

""^^~2(^+2»o-/-^-)' (37) 

(-1) _ -(no2i?2t {^a^ + n^RH-')) _ {2 + n^H^) 



no2 i;2 1 (lo - 1°"";,^^ + ^2 



= 4 + ^ (39) 

and 



(-3), , _ -{no'^ R^t {176 a^-l&Aa'^no'^ R^t^+i&O a? no'^ R^f^-2Sno^ R'^t'^)) 

^TM (*) = z r"9 1" 

128a«(l+ "°'f ) 

^ " ^ (40) 

_^t^ (2+no^ t^) (4+12 no^ t^+Sno* t*) Z^+2 (l+no^ t^)^ (8+12 rao^ t^+3no* **) Z* 

16 15 (l+no2t2)5 ■ 
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